Statistical properties of the quantum anharmonic oscillator. 
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I . Abstract 

The random matrix ensembles (RME) of Hamiltonian matrices, e.g. Gaussian random 
matrix ensembles (GRME) and Ginibre random matrix ensembles (Ginibre RME), are ap- 
plicable to following quantum statistical systems: nuclear systems, molecular systems, con- 
densed phase systems, disordered systems, and two-dimensional electron systems (Wigner- 
Dyson electrostatic analogy). A family of quantum anharmonic oscillators is studied and 
the numerical investigation of their eigenenergies is presented. The statistical properties of 
the calculated eigenenergies are compared with the theoretical predictions inferred from the 
random matrix theory. Conclusions are derived. 

II. Quantum harmonic oscillator in = 1 spatial dimension. 
Firstly: Let us consider the Hilbert space: 

Vi = L2(R,C,da;), (1) 

of the complex valued wave functions \E' that are (modulus) square integrable on the set 
R of the real numbers, and C is the set of the complex numbers. The Hilbert space Vi is 
separable space, and its orthonormal basis S is a set of Hermite's functions (Fok's func- 
tions, eigenf unctions of the dimensionless Hamiltonian operator 7i of the quantum harmonic 
oscillator in D = 1 spatial dimension): 

^^(x) = NrHr{x)exp{-^x^),Nr = [y/^r\2'']^ , r G N, (2) 
where N is a set of natural numbers including zero, whereas: 

H^ix) = (-l)-exp(x2)-^exp(-x2), (3) 

dx^ 

is rth Hermite's polynomial. We assume from now that all the operators are dimensionless 
(nondimensional). The dimensionless quantum operator x of the xth component of the 
position vector (radius vector) acts on the Hermite's basis function as follows: 




X<I/r = \h^r-l + \ (4) 



i 



whereas the quantum operator of the xth component of hnear momentum vector is given 
in the basis by the following formula: 

Matrix elements xi^^ and (pa;);,^ of these operators equal correspondingly: 



and 



where 




is discrete Kroneckcr's delta (it is not continuous Dirac's delta). The quantum operator 
of the square of the xth. component of radius vector after acting on the basis function reads: 

x^-^r = \ Wr + (2r + 1)*^ + Vr + iVr + 2*^+2], (9) 

and its matrix element {x^\,t is equal: 

{x\,r = \[VV^Mi,r-2 + (2r + l)5i,r + V7+l^^7+25l,r+2i (10) 

whereas the quantum operator of the square of the xth component of the linear momentum 
vector has the following representation in the basis: 

pI^v = ^[-Vr- lV^*,-2 + (2r + 1)*, - Vr + lVr + 2*,+2], (H) 
and its matrix element (p^)i,r is equal: 

{pDhr = \[-V^^V^5i,r-2 + (2r + 1)5,,, - V^^/^5l^r+2\. (12) 

Therefore, tie dimensionless (nondimensional) quantum Hamiltonian operator Ti. reads: 

n = pl + x\ (13) 

acting on the basis function, it produces: 

n^r = (2r + 1)^„ (14) 
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and its matrix element Hi^r is equal: 



ni,r = (2r + l)5i,r, (15) 

where 

= 2r + l, (16) 

is the rth eigenenergy of 7Y. The eigenenergies are simply all odd natural numbers, and the 
quantum Hamiltonian is diagonal operator, and its matrix representation is diagonal cxo x oo 
matrix. Note, that if one introduces the notion of nearest neighbour energy spacing (NNS): 

Sr = er+i - e„ r = 1, • • • (TV - 1), TV > 2, (17) 

then for the quantum harmonic oscillator it holds: 

Sr — 2 — const, (18) 

so the eigenenergies are equidistant. Here N is the number of energies dealt with (it is 
also a dimension of truncated Hilbert space Vi^n = C^). The eigenfunctions of Vi,Ar are 
A'^-component constant complex vectors (analogs of spinors), and the operators acting on it 
are N x N deterministic complex-valued matrices. The probability distribution Pjv-i of the 
spacings is discrete one point distribution for any finite value of N, N > 2: 

tending in the thermodynamical limit — > oo to the singular Dirac's delta distribution: 

Poo{s) = 6{s-2). (20) 

Secondly, let us perform more difficult task consisting of calculating all the absolute 
moments: _ 

{mn)l,r = mn= (x"),,, = {^l\x''^r) = / ^tix)x''^r{x)dx, (21) 

of the position operator x in the Hermite's (Fok's) basis B. One can calculate the lower 
moments manually, e. g., using recurrence relations, matrix algebra, etc., but it is tedious 
(even for 3 < n < 6). If one wants to calculate all the moments then he must return to the 
beautiful XIX century mathematics methods and after some reasoning he obtains the exact 
formula: 

{mn)l,r = "In = (a:'')«,r = 

[r/2] ^ 



^ ^ ^ ^t^o^t^o^ j!(/-2j)!^!(r-2^)! 
.2i+S-2,-2.-ir(" + ^ + r-^2j--2/c + l^j^ ^22) 
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where [•] is entier (step) function, F is Euler's gamma function. Therefore, the matrix 
representation of the even power operators x^'^ in the basis B are hermitean (symmetrical real) 
matrices with nonzero diagonal and nonzero p subdiagonals (and nonzero p superdiagonals), 
where the distance of the nearest superdiagonals (or subdiagonals) is 2 (the diagonal is also 
distant by 2 from the nearest super- and sub-diagonal), whereas the odd power operators 
^2p+i -j^ ^jjg basis B are hermitean (symmetrical real) matrices with zero diagonal and p 
nonzero subdiagonals (and p nonzero superdiagonals), where the distance of the nearest 
superdiagonals (or subdiagonals) is 2 (the nearest super- and sub-diagonal are also distant 
by 2). The physical interpretation of the superdiagonals (and subdiagonals) is connected 
with the absorption (emission) of phonons. 

III. Quantum anharmonic oscillator i.n D — 1 spatial 
dimension. 

Thirdly, we are ready to deal with the quantum anharmonic oscillator m. D = 1 spatial 
dimension. Its dimensionless Hamiltonian operator Tif^anharm reads: 

s 

'^l.anharm = 'H + X/ ^s^^ -i (23) 
s=0 

where is a degree of the anharmonicity of the oscillator, and the prefactors are the 
strengths of anharmonicity. The matrix elements of the anharmonic Hamiltonian operator 
are: 

5 

(^f,anharm)/,r = 6,5,,, + a.(^')^r-, (24) 

where the representation of the quantum anharmonic oscillator in the quantum harmonic 
oscillator basis B is mathematically correct, because the basis i3 is a complete set, and the 
Hilbert space of the eigenfunctions of the anharmonic oscillator is isomorphic to the Hilbert 
space Vi for the harmonic oscillator, provided that the anharmonic potential: 

s 

^l,anharm(^) — X/ ^s-^^' (^^) 
s=0 

is bounded from below (there are no scattering eigenstates) . 

Fourthly, we repeat the "Bohigas conjecture" that the fluctuations of the spectra of the 
quantum systems that correspond clasically to the chaotic systems generally obey the spectra 
of the Gaussian random matrix ensembles. The quantum integrable systems correspond to 
the classical integrable systems in the semiclassical hmit [1, 2]. The probabihty distributions 
Pp of the nearest neighbour spacing for the Gaussian orthogonal ensemble G0E(2) of 2 x 
2 Gaussian distributed real- valued symmetric random matrix variables {j3 = 1), for the 
Gaussian unitary ensemble GUE(2) of 2 x 2 Gaussian distributed complex-valued hermitean 
random matrix variables (/5 = 2), for the Gaussian symplectic ensemble GSE(2) of 2 x 2 
Gaussian distributed quaternion- valued selfdual hermitean random matrix variables {(3 = 4), 
and for the Poisson ensemble (PE) of the random diagonal matrices with homogeneously 
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distributed eigenvalues on the real axis R are given by the formulae: 

Pp{s) = e{s)Af,s^exp{-Bps'), (26) 
for the Gaussian ensembles, and 

Po{s) = e{s)exp{-s), (27) 

for the Poisson ensemble, where 6 is Heaviside's unit step function [3, 4, 5, 6, 7, 8, 9, 10, 11, 
12]. The constant numbers: 

_ r/^+H(/3 + 2)/2) _ r^((/3 + 2)/2) 

^'-^"mWTm) ^^-r2((/3 + i)/2)' (2^) 

are given by the formulae: Ai = 7r/2, Bi = tt/A (GOE), A2 = = 4/7r (GUE), 

and A4 = 262144/7297r'^, i?4 = 64/97r (GSE), respectively. For the Gaussian ensembles the 
energies are characterized by the "level repulsion" of degree (3 near the origin (vanishing 
spacing s = 0), and the probability distributions Pp vanish at the origin, and the quantum 
system with "level repulsion" is cast to the class of quantum chaotic systems. For the Poisson 
ensembles the energies are characterized by the "level clustering" near the origin (vanishing 
spacing s = 0), and the probability distributions Pq has maximum at the origin, and the 
quantum system with "level clustering" are treated as the quantum integrable system. After 
many numerical experiments conducted with different quantum anharmonic oscillators (up to 
the sextic quantum anharmonic oscillator 5" = 6) we draw conclusion that majority of them 
behaves like quantum integrable systems, the eigenenergies tend to cluster, the histogram of 
nearest neighbour spacing is closer to the Pq distribution resulting from the Poisson ensemble. 

IV. Quantum harmonic oscillator in D >1 spatial dimensions. 

Fifthly: Quantum harmonic oscillator in D spatial dimensions is a solvable analytically 
model. In order to make the deliberations easier we reduce our present interest to the first 
quantization. Therefore the relevant Hilbert space Vd is a D-dimensional tensor (Cartesian) 
product of the one-dimensional Hilbert spaces Vi: 

D 

Vi, = 0Vi, (29) 

whereas the harmonic oscillator's eigenfunction \E'r in D dimensions is, neither symmetrized 
nor antisymmetrized, tensor product of the eigenfunctions in one dimension: 

D 

*r(x) = n *r,(2^i),x = {xi, ...,xd) G R^, r = (n, ...,rD) e N^, (30) 
where we used boldface font for the £)-tuples x, and r. It follows that: 

D D 

^'r(x) = NrH^i^) exp(--x2), iVr = n ^r,, H,{^) = [] Hr.{Xj). (31) 
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One can also draw a conclusion that the Hilbert space: 

Vd = L\R^,C,dx),D>l, (32) 

is composed of the complex valued wave functions ^ that are (modulus) square integrable on 
the set R"^. The Hilbert space Vd is separable space, and its orthonormal basis Bd is a set of 
Hermite's functions in D dimensions (Fok's functions in D dimensions, eigenfunctions of 
the dimensionless Hamiltonian operator Hd oi the quantum harmonic oscillator in D spatial 
dimensions) . We assume from now that all the operators are dimensionless (nondimensional) . 
The dimensionless quantum operator Xj of the j'th component of the position operator x acts 
on the Hermite's basis function as follows: 



= + ^^-^^r,+l,j = 1,...,D, (33) 

whereas the quantum operator pj of the jth component of linear momentum operator p is 
given in the basis by the following formula: 



P^.^r, = - 7V^^^.+1- (34) 

Matrix elements {xj)i-^rj and {pxj)ij,rj of these operators are given by: 



and 



where 



^r=f[Sl^,r„ (37) 

is discrete Kronecker's delta in D dimensions (it is not continuous Dirac's delta). The 

of the square of the operator Xj 



quantum operator x^ of the square of the operator Xj after acting on the basis function 



reads: 



^Pr, = -[.Jr, - l,/¥]^r,-2 + (2r, + 1)^',^ + ^rj + l^rj + 2^r,+2], (38) 



and its matrix element {x^)ij,rj reads: 



1 



i^j)h,r, = 2 tV^^- - ^V^^h,r,-2 + (2r, + 1)5;^,,^ + ^Tj + l^Tj + 2Si^,r,+2]: (39) 

whereas the quantum operator p^. of the square of the operator p^. has the following repre- 
sentation in the basis: 

1 



Px.^r, = -[-^JrJ - lV^^'r,-2 + (2r,- + 1)^,, - ^r, + l^rj + 2^,.+2], (40) 
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and its matrix element ipl^)ij,rj is equal: 



= l[-^fr^V^^h,r,-2 + (2r,- + l)5i.,r, - ^Jrj + l^rj + 25i.,r,+2]. (41) 

Therefore, the dimensionless (nondimensional) quantum Hamiltonian operator TYd of the 
quantum harmonic oscillator in D dimensions reads: 

D 

nj, = p^ + y,^ = j2nj, (42) 

acting on the basis function, it produces: 

TY^^r = Tr(2r+l)*r, (43) 
and its matrix element {Tinj^r is equal: 

(^D)i,r = Tr(2r+l)5i,r, (44) 

whereas 

D 

er = Tr(2r+l) = ^e,,, (45) 
j=i 

is the rth eigenenergy of TYd. The eigenenergies are simply the sums of all odd natural 
numbers, and the quantum Hamiltonian is tensor product of diagonal operators in one 
dimension, and its matrix representation is tensor product of diagonal oo x oo matrices (it 
is poly-index matrix). There is no unique straightforward analog of the nearest neighbour 
spacing in D dimensions. 

Sixthly, let us perform very difficult task consisting of calculating all the absolute mo- 
ments: 

^ roo 

(m„)l,r = mn = (X")i,r = (^l|x"^r) = JJ / %{Xj)x]^-^rA^j)dXj, (46) 

of the position operator x in the Hermite's (Fok's) basis Bd- Here, as usual: 

i^ni,r-m^?h,r, (47) 

One can calculate the lower moments manually, e. g., using recurrence relations, matrix 
algebra, etc., but it is tedious (even for 3 < rij < 6). The exact formula for all the moments 
reads: 

D 

(mn)l,r = m„ = (x")i,r = Y[{mnj)lj,rj ■ (48) 
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V. Quantum anharmonic oscillator in D > 1 spatial dimensions. 

Seventhly, wc are ready to investigate the quantum anharmonic oscillator in D spatial 
dimensions. Its dimensionless Hamiltonian operator 7i|)_anharm reads: 

S {Si,...,Sd) D 

^Aanha^m = + E = + E ["{.i,..,.^.) " IK^)'']' (49) 

s=0 (si,...,si3)=(0,...,0) j=l 

where S is a D-tuple of degrees of the anharmonicity of the oscillator, and the prefactors 
are the strengths of anharmonicity. The matrix elements of the anharmonic Hamiltonian 
operator are: 

s 

(^Aanha.m)l,r = Cr^l.r + E «s(x^)l,r, (50) 

s=0 

where the representation of the D-dimensional quantum anharmonic oscillator in the quan- 
tum harmonic oscillator basis Bd is mathematically correct, because the basis Bd is a com- 
plete set, and the Hilbert space of the eigenfunctions of the anharmonic oscillator is iso- 
morphic to the Hilbert space Vd for the harmonic oscillator, provided that the anharmonic 
potential in D dimensions: 

S (5i,...,5d) d 

WB,anharm(x) = E «sX^ = E Ku-.so) ' (51) 

s=0 (si,...,sc)=(0,...,0) j=l 

is bounded from below (there are no scattering eigenstates in D dimensions). 

Eighthly, we repeat that the "Bohigas conjecture" also holds for the quantum oscillators 
in D dimensions. Having conducted many numerical experiments with different quantum 
anharmonic oscillators (up to the sextic {D = 3)-dimensional quantum anharmonic oscilla- 
tors: Sj = 6) we draw conclusion that some of them behave like quantum integrable systems, 
the eigenenergies tend to cluster, the histograms of nearest neighbour spacing are closer to 
the Pq distribution resulting from the Poisson ensemble, whereas other ones look like quan- 
tum chaotic systems, their eigenenergies are subject to repulsion, the histograms of NNS are 
closer to the distributions Pi, P2, P4, derived from the Gaussian Random Matrix ensembles. 
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